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field  on  the  stimulated  Brlllouin 


and  Raman  scattering  processes  are  studied.  The  formalism 
applies  to  plasmas  produced  by  CO^  lasers  and  to  electron 
cyclotron  heating  of  toroidal  systems  by  an  extraordinary 
electromagnetic  wave.  In  the  case  of  laser  fusion  the  plasma 
is  magnetized  due  to  the  self-generated  dc  magnetic  field 
while  in  toroidal  plasmas  it  is  due  to  the  external  magnetic 
field.  The  magnetic  field  greatly  reduces  the  threshold  for 


Brillouin  backscattering  by  the  lower  hybrid  wave.  The  Raman 


scattering  by  the  upper  hybrid  wave  has  substantial  growth 


kAn  because  of  the  lack  of  Landau  damping 
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I . INTRODUCTION 


The  absorption  of  the  laser  beam  by  plasma  is  dependent  upon  the 

scattering  processes  in  the  underdense  region.'*'  ^ As  resonance 

absorption  is  one  of  the  major  processes  by  which  the  laser  energy 

11  12 

is  deposited  in  the  plasma,  * the  laser  radiation  must  reach  the 
critical  density  layer  to  be  effectively  absorbed.  Hence,  it  is  necessary 
to  know  the  collective  phenomena  which  will  enhance  the  scattering 
of  the  incident  radiation  and  prevent  the  electromagnetic  wave  from 
reaching  the  critical  density. 

13-15 

It  has  been  measured  experimentally  that  in  laser-pellet 

experiments  there  is  a self-generated  dc  magnetic  field  with  its 

direction  perpendicular  to  the  polarization  and  propagation  vectors 

of  the  incident  laser  radiation.  Magnetic  field  of  several  megagauss 

16—22 

has  also  been  predicted  theoretically  and  observed  in  computer 

23 

simulations.  The  magnetic  field  is  generated  at  the  critical  layer 
and  spreads  out  to  the  less  dense  regions  of  the  plasma  by  convection 
and  diffusion.  Therefore,  the  investigation  of  the  stimulated  scattering 
processes  in  a magnetized  plasma  is  a relevant  part  in  the  theoretical 
understanding  of  the  fusion  driver-matter  interactions. 

24 

With  the  development  of  high-power  cyclotron  masers  the  electron 
25-30 

cyclotron  heating  becomes  an  Important  plasma  heating  mechanism 

31  32  33 

of  toroidal  systems  like  tokamaks,  ’ Elmo  bumpy  torus,  etc.  If 
the  extraordinary  electromagnetic  wave  is  launched  into  the  plasma 
from  the  side  of  the  lower  magnetic  field  it  needs  to  tunnel  through 
the  evanescent  region  produced  by  the  cyclotron  cutoff  in  order  to  be 
linearly  converted  into  the  electron  Bernstein  wave  at  the  upper 
hybrid  resonance.  In  addition,  the  electromagnetic  energy  can  be 
transferred  to  the  plasma  electrons  by  cyclotron  damping  at  the 


1 o.v 
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electron  cyclotron  resonance.  However,  if  the  extraordinary  electro- 
magnetic wave  is  injected  into  the  plasma  in  the  direction  of  increasing 
major  radius  the  cyclotron  cutoff  region  is  avoided. 

Similar  to  the  problem  found  in  laser  fusion,  the  stimulated 
scattering  of  the  incident  electromagnetic  wave  may  prevent  the  energy 
from  reaching  the  upper  hybrid  and  electron  cyclotron  resonance  layers. 
Therefore,  in  the  present  paper  we  also  examine  the  role  played  by  the 
stimulated  Brillouin  and  Raman  scattering  instabilities  in  reflecting 
the  incident  microwave  energy  launched  into  the  toroidal  plasmas  for 
the  purpose  of  supplementing  the  initial  Joule  heating. 

For  the  stimulated  Brillouin  scattering  in  a magnetized  plasma, 
the  enhanced  scattering  is  due  to  the  excitation  of  an  electrostatic 
lower  hybrid  wave  instead  of  the  usual  ion  acoustic  wave  in  an  unmag- 
netized plasma,  while,  in  the  stimulated  Raman  scattering  it  is  due  to 
the  upper  hybrid  wave  instead  of  Langmuir  waves. 

In  Sec.  II,  the  general  formalism  of  nonlinear  scattering  pro- 
cesses in  a magnetized  plasma  is  developed.  In  Sec.  Ill,  the  growth 
rates  for  the  stimulated  Brillouin  and  Raman  scattering  are  found  with 
applications  to  fusion  produced  by  COj  lasers  and  electron  cyclotron 
heating  in  toroidal  devices.  Thresholds  and  the  effects  of  collisionless 
damping  on  the  Raman  instability  are  discussed  in  Sec.  IV.  In  Sec.  V, 
the  discussion  and  conclusions  are  presented. 

II.  GENERAL  FORMALISM 

Consider  an  extraordinary  electromagnetic  pump  wave  ^“(E^ad-Eg^y)* 
cos(kgX-Wgt)  incident  on  a plasma  with  the  static  magnetic  field  Bq 
along  the  z direction.  The  Lorentz  force  drives  the  electrons  with 


density  oscillations 
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1 

l 


and  e,m  being  the  electron  charge  and  mass,  respectively,  c is  the  speed 

i 

of  light. 

The  purpose  of  the  present  work  is  to  study  the  interaction  of  this 
pump  wave  with  a scattered  electromagnetic  wave  and  an  electrostatic 
decay  waVe,  whose  Fourier  components  are  (k_*k-kQ,  io_**u)-Wq)  and  (k,b>), 
respectively.  The  electromagnetic  waves  are  subjected  only  to  a 
collisional  damping  which  is  weak  and  can  be  neglected.  Although  fluid 
equations  will  be  used  in  the  present  formalism,  the  collisionless 
damping  of  the  waves  can  be  taken  into  account  phenomenologically. 


5 


The  equation  for  the  electrostatic  wave  is  determined  by  using 


Poisson's  equation 


v $*4ne(n-nj) 


where  <f>  is  the  perturbed  electrostatic  potential  and,  n and  n^  are 
the  perturbed  electron  and  ion  densities,  respectively.  The  wave 
equation  for  the  electromagnetic  decay  wave  is  found  from  Maxwell's 


equation  as 


2 ^ 

V2E_-V ($•£_)  + — ~ eL  •£_ 


where  ^ is  the  perturbed  electric  field  and  3^  is  the  nonlinear 
part  of  the  perturbed  current  density  originating  from  the  heating 
of  the  oscillating  velocity  in  the  incident  electromagnetic  wave 
with  the  low  frequency  density  oscillations 

^*?L-en?_  " \ en^  - 1 enK  * (8) 

The  anti-Stokes  component  was  neglected  as  being  off-resonant. 

The  dielectric  tensor  is  defined  as 

| =f+4nij  (9) 

L_  o)_  - 

where  f is  the  unity  tensor  and  o is  the  conductivity  tensor  related 
to  the  linear  current  density  by 

jk*o  . (10) 

The  perturbed  density  oscillation  n and  velocities  V and  are  determined 
from  the  continuity  equation  and  the  equations  of  motion,  respectively: 

|a  + ?•  (njj»  + 1 «0?_  + \ n.V-o  • (11) 


6 
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3V_  1 ->*  -fc*.  1 -►  ±+*  eE- 

ir  + 2 Vo'^  + 2 ^’^0  ■ ” “m“ 

4x3  -^-Vx3* 

- e 2mc  0 


where  V2*T/m,  T being  the  electron  thermal  energy,  and  B_  is  the  perturbed 
magnetic  field  of  the  scattered  wave.  The  perturbed  electron  density 
can  be  separated  into  the  linear  and  nonlinear  parts 


L.  NL 
n-n  +n 


where 


!V  (k2  + _»i_  k2\ 

» ("  „2-n2  V 

e 

2 (.2  . u.2  .2  \ 2 

“ ‘ lkll  2 02  iVe 

\ to  -n  ' 


k(|  and  kx  designating  the  parallel  and  perpendicular  components  of 
k in  relation  to  the  static  magnetic  field,  and  n^1  is  found  from  Eqs. 
(11) — <13) . The  ions  are  considered  cold  and  unmagnetized.  The  nonlinear 
effects  such  as  the  ponderomotive  force  on  the  ions  can  be  neglected 
as  it  is  smaller  by  a mass  ratio  to  that  of  electrons.  The  ion 


perturbed  density  is  given  by 


0, 2. 

enQk  + 


where  M is  the  ion  mass.  The  substitution  of  Eqs.  (14)-(16)  into  Eq. 
(6)  gives 


(17) 
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where  eT  is  the  dispersion  function 


2 „2ft-  + -f 

L 2 o / o 

to  2 / i,2  . to 

to  - I + -r 


/ki.  . „2  k*\ 
-(kn  + -5^k‘) 


Equation  (15)  can  be  rewritten  as 

eT  k(t>*ot^*E 

where  is  expressed  by  the  transpose  of  matrix  a describing  the  low 
frequency  beating  of  the  two  electromagnetic  waves.  After  Fourier 


analyzing  Eq.  (7),  we  get 


=f-(tc<f>) 


where  1$  is  the  usual  dispersion  tensor 

?-k2i+££ 


and  B is  expressed  by  the  matrix  containing  the  coupling  between  the 
pump  wave  and  the  low  frequency  electrostatic  wave.  The  components 
of  the  dielectric  tensor  are  found  from  Eqs.  (9)-(13)  to  be 


2 2 
io  -ft 

- e 
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ft  (0 

-i  _i— E— 
u)  2 2 
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Equations  (19)  and  (20)  comprise  the  set  of  coupled  equations  describing 
the  parametric  processes  under  investigation. 

In  the  following  we  shall  discuss  the  scattering  in  the  plane 
perpendicular  to  the  magnetic  field  which  is  also  the  maximum  growing 

modes  for  backscattering.  In  this  case  k|(  =0  and  Eq.  (18)  reduces  to 

2 2 
m . u> 

el“1  2 2 — 2?  2 2 (23) 

L u>  u-Sl-kV 
e e 

, -+T 

and  the  nonvanishing  components  of  a are 


and 


2o2 

b)  0 

JLJL 


2 2 2 2 2 2 

2(B-f)(/-i)-kT) 
- e e e 
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mfi  mm-mfl 

u e 
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- e e e 
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e 
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0) 
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(24) 
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From  Eq.  (21)  the  components  of  the  dispersion  tensor  become 


(25) 


a)2  n 


(U 


2 o2 
u)  -0  , 
- e 


m 


H _ £ — 

2 m 2 _2 
c - m -0 
- e 


2 „ 2 
m (2  m 
i — -£  — P- 
2 w 2.2 
c - m -ii 
- e 


2 / 2 \ 
m / o>_  \ ~ 

“if1"  2?2  )~  k- 


2 , 

m_  / m 

-i(i-- 

c \ m 


2 


- k‘ 


(26) 
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and  the  relevant  components  of  the  coupling  tensor  6 are 


* 

ik  V 
0 Ox 


„ 2,  2 2W  2 2 ,2  2, 
2c  (u  -0  ) (w  -0  -k  V ) 
- e e e 


0)qO)_  2 2 
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Expressions  (24) -(28)  reduce  to  those  of  unmagnetized  plasma  If  we 


set  0 =0. 
e 


From  Eq.  (20)  we  obtain 


Ad^1 


t.-tr1*  [?•  <£♦)  i = -f]-^n — • [$•(&!>)] 


where  ||d||  represents  the  determinant  of  the  matrix  D and  Ad^  Ir1 
denotes  the  adjunt  of  the  inverse  to  the  matrix  D.  Substituting 
Eq.  (29)  into  Eq.  (19)  we  get  the  expression  for  the  dispersion 
relation 


HdJI  = aT*  [AdjS^1*  ($•£)  ] 


where 


10 


t-1  _ i-1 


D D -D  D 
xx  yy  xy  y 

D 


, (31) 


D.  = D D -D  D , (32) 

11  1"  xx  yy  xy  yx 

and  k=K/k.  The  dispersion  relation  (30)  is  valid  for  the  scattered 

electromagnetic  wave  making  an  arbitrary  angle  with  the  incident  wave. 

It  also  applies  if  the  electrostatic  wave  is  either  an  eigenmode  or  a 

quasimode.  If  we  assume  resonant  scattering  by  eigenmode,  i.e.,  e^(k,m)=0, 

we  can  expand  e and  ||d  ||  about  the  eigenf requencies  u>  [e  (w  )=0]  and 
L 1 k L k 

0)  ^ [ || D ||  (a)  ^ )=0]  to  get  the  growth  rate 

~k_  ~k_ 

2 aT* [AdjD^1* (f *k) ] 


For  the  backscattering  process,  expression  (33)  reduces  to 


0 D a 6 +D  a 6 
2 yy  x xx  xx  y i 


Y=_8eL 
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For  a scattering  process  in  the  plane  perpendicular  to  the  magnetic 
field  k(|  »0,  JN^  vanishes  and,  hence,  0zx=O.  Consequently,  the 
decay  of  an  extraordinary  electromagnetic  wave  into  an  electrostatic 


and  an  ordinary  electromagnetic  waves  is  not  possible. 

For  a nonresonant  scattering  by  the  quasimode,  the  growth  rate  is 


given  by 


iaT* [AdjD^1* (f*k) ] 


EL  8w 


For  the  backscattering  case  this  expression  reduces  to 


Y = - 


i(D  a 6 +D  a 0 ) 

yy  x xx  xx  y yx 
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III.  BRILLOUIN  AND  RAMAN  SCATTERING  INSTABILITIES 


A.  Laser  Fusion 

The  validity  of  this  investigation  is  limited  by  the  fluid 

assumption  which  for  a magnetized  plasma  requires  that  ^k^p^O.S  (pg= 

3 -1 

V /ft  ) . For  instance,  for  a C0o  laser  (k-.=5.9xl0  cm  ),  a self- 
e e l U 

13 

generated  magnetic  field  of  2 Mgauss  (fie=3.5xl0  rad/s),  and  plasma 

9 

temperature  of  1 keV  (Ve=1.3xl0  cm/s),  then,  2kgPe=0.4  and  the  theory 
is  only  valid  marginally  for  these  parameters.  Obviously,  for  CO2 
lasers  and  stronger  magnetic  fields  or  colder  plasmas  the  condition  of 
validity  of  the  theory  is  well  satisfied.  In  the  case  of  Nd-glass 
laser  (^=5.9*10  cm  ),  for  a 4 Mgauss  magnetic  field  (0e=7xl0  rad/s), 


the  theory  is  valid  only  in  a very  cold  plasma  with  maximum  plasma 


12 


O 

temperature  of  50  eV  (Vg=3xl0  cm/s).  We  conclude  that  kinetic  theory 
should  be  used  to  describe,  more  realistically,  plasmas  produced  by 
Nd-glass  lasers  where  the  temperatures  are  the  order  of  a few  keV. 

In  the  kinetic  formalism,  it  is  possible  to  study  the  scattering  by 
Bernstein  modes  which  might  give  quite  different  results.  This 
remains  to  be  investigated. 

To  study  the  stimulated  Brillouin  backscattering  in  a magnetized 

plasma  in  the  limit  valid  for  laser  fusion,  a)p>>fie,  the  low  frequency 

1/2 

normal  mode  is  obtained  from  Eq.  (23)  as  w-(n  ft.)  • By  reducing 

ei 

the  expressions  (24),  (25),  (27),  (28),  and  ( 35) — ( 37 ) to  this  limit 

and  substituting  into  Eq.  (34),  the  growth  rate  becomes 

1/2 

i ...  I Pi  I / n 

YB 


-Wofe)  (if 


(40) 


The  growth  rate  in  an  unmagnetized  plasma  is  given  by 

o 1/2 


U J 

_eL 


— k V 
2 0 0 Wkcs 


(41) 


The  ratio  between  these  two  growth  rates  is  given  by 


f " (kpe) 


1/2 


(42) 
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For  instance,  if  Wq=1. 78x10  rad/s  (CO2  laser),  Bg=2  Mgauss,  and 

T =1  keV,  then  y /y-O.bb. 
e d 

For  the  stimulated  Raman  instability,  Eq.  (23)  reduces  to 
2 2 1/2 

(D*(u  +ft  ) . From  Eq.  (34)  the  growth  rate  is  found  to  be 

p 6 ,1/2 

1 / \ 

v s — k V 

In  o i'tv  v , 


2 0 0 Vv>uh 


(43) 


The  growth  rate  for  the  Raman  scattering  in  an  unmagnetized  plasma  is 

.1/2 

Y - -y  knv„  | — I . (44) 


2 Vo  \»0; 


T 


which  is  close  to  one  for  laser  fusion  parameters. 

B.  Electron  Cyclotron  Heating 
The  condition  2kQpg<0.5  for  the  validity  of  the  present 
formalism  is  well  verified  for  microwave  injection  in  toroidal  plasmas 
in  the  direction  of  increasing  major  radius.  The  magnitude  of  kg 
increases  from  zero  at  the  upper  hybrid  cutoff  to  the  vacuum  value 
Wq/c  at  the  plasma  boundary.  Supposedly,  we  want  to  deposit  the 
wave  energy  in  the  central  part  of  the  system,  then,  the  upper  hybrid 
resonance  should  be  as  close  as  possible  to  the  center  of  the  plasma 
and  that  is  achieved  by  choosing 

a)  =[a,2(r=0)+ft2(r=0)  ]1/2  (46) 

0 uh  p e 

where  r is  the  variable  which  describes  the  minor  radius  of  the 
toroidal  device.  If  we  assume  Wp(r=0) /Qe(r=0)=!l , then  Eq.  (46) 
gives  o)Q=»,/2ne(r=0) . The  toroidal  magnetic  field  in  tokamaks  decreases 
as  1/R  when  R increases  where  R describes  the  major  radius,  then,  the 
maximum  value  of  the  electron  cyclotron  frequency  in  the  region  where 
the  electromagnetic  wave  propagates  occurs  just  before  the  upper  hybrid 
cutoff  layer.  If  we  assume  a parabolic  density  profile  to  determine 
the  cutoff  layer,  we  find  that  fie(r«cutoff)=0.9  (r=0)  and,  then. 


p =V  / 0.9  n (r=0) . We  conclude  that 

KP  P P 


— 
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For  microwave  injection  from  the  inner  part  of  the  torus,  the  theory 

is  not  valid  for  Brillouin  and  Raman  scattering  processes  near  the  upper 

hybrid  resonance  layer  since  the  refraction  index  becomes  infinity 
2 2 1/2 

when  u>  =(u>  +fi  ) . However,  in  practice  we  like  also  to  have  the 

1)  p e 

wave  energy  transferred  to  resonant  electrons  by  cyclotron  damping 

if  we  inject  the  microwave  with  finite  k(|  . In  this  case,  u)Q“fie+k()  V(( 

and  the  electromagnetic  wave  reaches  the  cyclotron  resonance  layer 

before  reaching  the  upper  hybrid  resonance.  Then,  for  , 

k„=wri/c=ft  / c,  and  we  have 
u u e 


2V 


a v 

2k_p  =2  — ^ « — <<  1. 
0 e c Si  c 
e 


(48) 


To  investigate  the  Brillouin  backscattering  process  in  the  limit 

Wp=Qe,  valid  for  electron  cyclotron  heating  of  toroidal  plasmas, 

the  frequency  of  the  electrostatic  wave  obtained  from  Eq.  (23) 

becomes  u)=u  . . The  growth  rate  is  obtained  by  reducing  the  expressions 
pi 

(24),  (25),  (27),  (28),  and  (35)-(37)  to  this  limit  which  by 
substitution  into  Eq.  (34)  gives 


'B 


2 , 1/2 
I k V i 


(49) 


The  ratio  between  the  growth  rates  applicable  to  laser  fusion 


and  electron  cyclotron  heating  is  given  by  the  ratio  of  the  expressions 

3 

(40)  and  (49),  respectively,  which  is  . 

For  the  stimulated  Raman  scattering  instability  the  growth  rate  is 
the  same  as  that  given  by  Eq.  (43).  The  ratio  between  this  growth 
rate  and  the  one  for  an  unmagnetized  plasma  is  given  by  Eq.  (45) 


which  is  2 


,-iM 


for  toroidal  plasma  parameters. 


15 


IV.  DAMPING  AND  THRESHOLDS 


The  collisionless  damping  for  upper  hybrid  waves  Is  given  by 

2 . 


2 2 2/ 

1/2  (AkV 

it  1 p _e 

2 2 2 k V 

k Xp  KH  e 


0), 


. 1 .2  2 

+ 2 ktPe 


exp 


(- 


(lvV 

k2  V2 
li  e 


exp 

2 > 


2 2 
k7.  V* 
II  e 
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(50) 


where  w2=co2+«2+k2V2.  For  scattering  instability  in  the  plane  perpendicular 
k p e e 

to  the  magnetic  field,  k,,  vanishes  and,  consequently,  the  damping  is 

4 

greatly  reduced.  In  the  case  of  unmagnetized  plasmas,  however, 

the  Langmuir  waves  become  heavily  Landau  damped  for  kX^>0.2  (l^V^/w^) 

and  the  Raman  scattering  instability  is  turned  off,  although,  the 

nonlinear  Compton  scattering  of  the  heat  mode  still  persists  with  a 

smaller  growth  rate.  But  for  magnetized  plasmas,  Raman  mode  prevails 

over  nonlinear  Compton  scattering  even  for  large 

If  we  assume  a weakly  inhomogeneous  plasma,  the  Rosenbluth 
34 

threshold  condition 


V1V2K' 


> 1 


(51) 


should  hold  where  Yq  is  the  growth  rate  in  a homogeneous  plasma, 
and  are  the  components  of  the  group  velocities  of  the  decay  waves 
along  the  density  gradient,  and  K * ™d/dx^ ^^oX— k2x^  C^e  der*vat*ve 
in  relation  to  the  spatial  variable  along  the  direction  of  the  density 
gradient.  For  the  Brillouin  scattering  instability,  the  group  velocity 
of  the  lower  hybrid  wave  is  given  by 


Hh 


(52) 
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Therefore,  for  scattering  process  in  the  plane  perpendicular  to  the 
magnetic  field  the  group  velocity  vanishes  and  from  Eq.  (51)  we  conclude 
that  there  is  no  inhomogeneity  threshold  for  this  process.  For  the 
Raman  scattering  process,  the  group  velocity  of  an  upper  hybrid  wave 
is  given  by 


g . ' 2u>  k'L 
°uh  uh 

which  substituted  into  Eq.  (51)  together  with  the  Raman  homogeneous 
growth  rate  for  magnetized  plasma,  gives  the  same  threshold  as  in  the 


unmagnetized  case 


i(?)v 


where  dK/dx»dk/dx  for  Raman  scattering. 


V.  DISCUSSION  AND  CONCLUSION 


In  the  fluid  theory  limit,  the  growth  rates  for  Brillouin 

scattering  process  in  magnetized  plasmas  are  reduced  in  relation  to  the 

growth  rates  for  unmagnetized  plasma.  For  plasmas  produced  by  CO^ 

1/2 

lasers,  the  growth  rate  is  reduced  by  a factor  (kpg)  while  for 

1/2 

electron  cyclotron  heating  of  toroidal  plasmas  by  a factor  (kpe)  * 

3 

(w  /0  ) . The  Raman  growth  rates  are  practically  unmodified  as  compared 
P ® 

to  the  unmagnetized  plasma  either  for  laser  fusion  or  electron 
cyclotron  heating  parameters. 

The  collisionless  damping  is  vanishingly  small  for  stimulated 
Raman  scattering  process  in  a plane  perpendicular  to  the  magnetic 
field  and,  consequently,  Raman  Instability  is  the  prevailing  process 
over  the  nonlinear  Compton  scattering  for  all  values  of  kX^. 
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I 
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By  assuming  a weak  inhomogeneity  in  the  plasma  density,  we  find 
that  there  is  no  threshold  for  the  stimulated  Brlllouin  scattering 
instability  in  a plane  perpendicular  to  the  magnetic  field.  For 
Raman  scattering  the  threshold  is  the  same  as  that  in  an  unmagnetized 
plasma. 
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